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ABSTRACT 

We report VIMOS integral field spectroscopy of the N159F nebula surrounding LMCX-1. Our 
observations reveal a rich, extended system of emission line filaments lining the boundary of a large 
conical cavity identified in Spitzer mid-IR imaging. We find that X-ray photoionization cannot be 
solely responsible for the observed ionization structure of N159F. We propose that the extended 
filamentary emission is produced primarily by ionization from a shock driven by a presently unobserved 
jet from LMCX-1. We infer a shock velocity of v s » 90kms~ 1 and conclude that the jet responsible 
for the bow shock is presently undetected because it has switched off, rather than because it has a 
low surface brightness. This interpretation is consistent with the present soft X-ray spectral state of 
LMCX-1 and suggests the jet is intermittent. 

Subject headings: accretion — ISM: jets and outflows — techniques: spectroscopic — X-rays: binaries 
-X-rays: individual (LMCX-1) 



1. INTRODUCTION 

LM CX-1 resides i n the highly complex nebula region 
N159 (Hcnizc 195(f). It is the most powerful known 
X-ray source in the LMC (jJohnston. Bradt fc Doxsevl 

with an X-ray luminosity Lx ~ 2 x 10 ergs 
( Schlegel et al.l I1994D . and is believed to be pow- 
ered by accretion onto a stellar mass (4 — 10 Mq) 
black hole wit h a companion 07ili-type star ident i- 
fied as star 32 (iPakull fc Angebauld fl^lCowlevlfl992l). 
iPakull fc Angebaultj (|1986l ) collated extensive spectral 
data and ruled out collisional excitation by a supernova 
remnant. They concluded that their detection of Hen 
emission is the first evidence f or an X-ray ionized nebula 
(see also iRamsev et all 120061 ) . We note, however, that 
the slit spectra used in these studies were limited to a 
relatively localized region of the nebula, in the immedi- 
ate vicinity of LMCX-1. 

Using integral field spectroscopy (IFS), we present new 
evidence suggesting that X-ray photoionization cannot 
be solely responsible for the observed large scale mor- 
phology and high excitation lines, which we propose are 
best understood in terms of shock ionization attributable 
to an as yet undetected relativistic jet from LMCX-1. 
Indeed, the highly complex ionization structure of the 
extended nebula surrounding LMCX-1 is strikingly rem- 
iniscent of the other two only known examples of jet- 
induced shock ionized nebulae, those su rrounding the 
powe rful Galactic X-ra y binaries SS433 (jDubner et al.l 
Il99"8h and Cygnus X-l (jOallo et alJl200fih . 

We were initially alerted to the possibility of shock 
ionization in N159F after observing with the Taurus 
Tunable Filter (TTF) at the Ang lo-Australian Telescope 
(|Bland-Hawthorn fc Joneill998f ) in 2003 December. We 
discovered a rich system of streaming Ha filaments in 
the arcmin-scale nebula. The filaments trail back from 
the apex of a conical cavity in the dust emission seen 
in the Spitzer Infrared Array Camera (IRAC) image of 
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iJones et al.l (|2005f ). The apex is positioned m 12" South- 
West of LMCX-1, correspon ding to ~ 3 .2 pc at a dis- 
tance d w 55kpc to the LMC (|Feastjll999f ). 

Motivated by these observations, we undertook further 
investigation into the environment of LMCX-1 through 
integral field spectroscopy (IFS) with VIMOS at the Very 
Large Telescope (VLT) in 2005 November. The new data 
reveal strong detections of optical emission lines due to 
[Oi], [Nil], [Sn], He I and [Arm], as well as Ha. The 
emission lines all show a similar morphological bow shock 
structure as the initial Ha filam entary streamers seen in 
the TTF image. In Figure 3 of I Jones et all (l2005lL which 
overla ys the 5 GHz radio contours of iHunt fc Whiteoakl 
(|1994h on an IR image of N159, radio emission appears 
to trace the IRAC cavity. 

We summarize the observations and data reduction 
procedures in § 2. In § 3, we present our results and pro- 
vide an interpretation of the environment of LMCX-1. 
Our conclusions are summarized in § 4. 

2. OBSERVATIONS AND DATA ANALYSIS 

Our initial TTF observations of the N159F nebula 
were taken in a 5 A bandpass centred on Ha. A 
powerful differential technique calle d 'straddle shuffling' 
(jMalonev fc Bland-Hawthorn1l2001h was used to remove 
the surrounding continuum emission to better than 1% 
accuracy, revealing the filamentary structure of the neb- 
ula. These data prompted observations with the VIMOS 
integral field unit (IFU) at the VLT on 2005 November 
2, 5, 6. Conditions were reported as photometric, with 
1.4" — 1.5" seeing (well matched to the 0.67" spatial scale 
of the IFU). Observations were queue scheduled using 
the HR-Orange grating (R«2500, AA5250-7400A) and 
include the emission lines Ha, [Oi], [Nil], [Sn], Hei, and 
[Anil]. The systemic radial velocity (v r ~ 270kms~ ) 
of the LMC ensures that the [O i] and Ha emission lines 
are well separated from their telluric counterparts. 

Given the 27" x 27" field-of-view of the VIMOS IFU 
in high resolution mode, four VIMOS pointings were 
required (with a 5 pixel/3.3" overlap) to tile the Ha 



2 



Cooke et al. 



filaments along the boundaries of the Mid-IR cavity 
in the vicinity of LMCX-1. Three 1160sec exposures 
were taken at each of the four pointings. The data 
were reduced using t he reduction pipeline VIPGI 3 (see 
iScodeggio et al.ll2005l for details). Datacubes were cre- 
ated and sub-sections of the final mosaic were combined 
using a set of IFU data manipulation routines created in 
IDL. Emission line fitting was performed for each promi- 
nent spectral line using a single unresolved Gaussian pro- 
file. 

Due to the restricted spectral coverage of the chosen 
VIMOS setting, the key shock diagnostic line [O m]A5007 
is not present in our data. However, we retrieved archival 
[Oni] observations from the public ESO archive. The 
2.3 m WFI observations 4 of N159F consists of 2 x 500 sec 
observations using the [Om] filter and 2x125 sec ob- 
servations using the R c filter. The data were pro- 
cessed in the usual manner (overscan corrected, flat- 
fielded from twilight flat frames, aligned and com- 
bined) using elements of the CASU imaging processing 
toolkit (jlrwin fe Lewis! [200 Continuum subtraction 
was achieved for the [O in] and Ha narrow band images 
using the aligned and scaled V and R c images respec- 
tively. The WFI and VIMOS data were spatially reg- 
istered and resampled using stellar images common to 
both data sets. 

3. RESULTS & DISCUSSION 

Figure 1 shows the main results from the new VI- 
MOS, WFI and IRAC observations. Together, these 
provide compelling evidence for shock ionization in the 
N159F nebula surrounding LMCX-1. Fig. l(a-d) show, 
respectively, the prominent emission lines Ha, [Ar in] and 
[O i] in the VIMOS data cube, and [O III] from the con- 
tinuum subtracted WFI data. A bow shock morphol- 
ogy, with trailing filaments, is evident in the images. 
Fig. 1(e) shows a line diagnostics plot formed from a 
pixel-pixel comparison of [0iii]/[0i] vs. [Ariii]/[Sn] 
where four color-coded ionization regions show spatial co- 
herence over the N159F nebula (see Fig. 1(f)). Fig. 1(h) 
is a 3-color Spitz er IRAC image centered on star 32 
([Jones et al.ll2005[ ) , showing a dust cavity associated with 
the proposed bow shock. Fig. l(k) presents a 3-color 
composite image in [N n] (red) , Ha (green) and [O i] 
(blue). We speculate that the shock is driven by a jet, 
with orientation indicated in Fig. l(k). 

3.1. Disentangling the Shock-Ionized Nebula 

Separating out the shock-ionized part of the opti- 
cal nebula is not straightforward, as radiative shocks 
can produce collisionally excited optical emission lines 
similar t o those produced by X-r ay photoionization 
(see e.g. iDopita fc Sutherland! [2003]) . Both ionization 
mechanisms can produce enhanced (relative to Hll 
regions) Ha as well as enhanced forbidden lines of 
neutral or low ionization species, such as [Oi]A6300, 
[Sn]AA6716,6731 and [Nn]AA6548,6583, and of higher 
ionization species, such as [O in] AA5007, 4959 and 
[Ann]A7136. HeiiA4686 is a reliable X-ray photoion- 
ization diagnostic as it is difficult to produce with 
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shock ionization, while [Ann] is difficult to produce 
with X-ray photoionization, except for very high ioniza- 
tion parameters (R. Sutherland, private communication) . 
Furthermore, X-ray photoionization tends to favour 
the p roduction of [Oni] over [Oi] ([Draine fc McKed 
1993). On the other hand, [Om] can be a power- 
ful shock diagnostic, s i nce all radiative shock models 
(ICox fc Ravmond|[i 985: Bi nette. Dopita fc Tuobyl 119851 ; 
lHartigan. Raymond fc Hartmannlll987j ) consistently pre- 
dict a dramatic onset of [O in] emission near a critical 
shock speed, v s ~ 100 km s" 1 . 

Whilst the Ha and [Arm] images in Fig. 1(a) and (b) 
exhibit a similar morphological structure suggestive of a 
bow shock, a comparison of the [O i] and [O in] images in 
Fig. 1(c) and (d) shows regions of anticorrelation in the 
imme diate vicinity of LMCX-1 and around the B5 1 star 
R148 (jFeast, Thackeray fc Wesserlinkll96"ol) . where there 
is a deficit of [O i] emission, but enhanced [O in] emission. 
R148 is not capable of producing this high excitation 
emission. Indeed, its Stromgren sphere (r s « 0.04 pc) 
does not even encompass one pixel of our image. The 
Stromgren radius for star 32 (the 07 in companion to 
LMCX-1 - ICui et alJ l2002T > is w 0.9 pc, comparable to 
the size of the blue region defined in Fig. 1(e) and (f). 
Recalling that X-ray photoionization favours [O ill] over 
[O i] and can produce [Ar in] for high ionization param- 
eters, this blue region, which exhibits the largest lo- 
cal enhancement in [Om]/[Oi] and [Ar in]/ [S n], must 
be due to the combined effects of X-ray ionization by 
LMCX-1 and photoionization by star 32. The yellow 
region in Fig. 1(e) and (f) is consistent with the extent 
of the X-ray Stromgren ra dius (~ 2pc), calculated by 
iPakull fc Angebaulti (| 19861 ). Note, however, that this 
yellow region is located to the East (left) of R148, on the 
far side of LMCX-1 and spatially coincides with some 
of the streamers. Thus, we interpret the yellow zone 
of enhanced [Arm]/[S n] and [Om]/[Oi] as being due 
to a combination of shock ionization and X-ray ioniza- 
tion. The red zone in Fig. 1(e) and (f) is predominantly 
shock-ionized. Fig. l(i) and (j) show a proposed trail of 
streamers through this shock-ionized region. 

The evidence for X-ray photoionization in N159F 
has hitherto been based on slit spectra of He n A4686 
e mission detected in the immediate vicinity of LM C X- 
1 (|Pakull fc Angebaulti [l98l iRamsev et alJ |2006[) . In 
Fig. 1(1) we show the posit ioning of the echelle slits 
used by Ramsey et al. (20Q 6[) to detect the Hell emis- 
sion. IPakull fc Angeba ult (198(f) used similarly placed 
orthogonal slits, but centred on star 32 instead of R148. 
Note that although the echelle slits coincide with the yel- 
low zone we identify as being both shock ionized and X- 
ray ionized, the observed He II A4686 emission is almost 
certainly due to only X-ray ionization as none of the ra- 
diative shock models predict significant He II A4686. Note 
also that it would have been extremely difficult to find 
unequivocal evidence for shock ionization from other line 
diagnostics due to the confusion with X-ray ionization in 
this region and also to the spatially limited extent of the 
slits. On scales of several parsecs, our IFU data reveal en- 
hanced [Oi]/Ha, [Om], and [Arm], as well as Ha, form- 
ing a shell-like structure around LMC X-l, but extending 
beyond the X-ray Stromgren radius (that is, beyond the 
yellow zone in Fig. 1(e) and (f)). This is compelling evi- 
dence that shock-ionization rather than X-ray ionization 
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is primarily responsible for this high-excitation emission 
seen on the largest scales in the nebula. 

3.2. Physical parameters 

Fig. 1(g) shows an image of the electron num- 
ber density, N e , deduced from the [Sn] emission 
(|Osterbrock fc FerTa nd 2006]). It reveals an enhanced 
flattened region that falls slightly inside and near the 
apex of the bow-shock morphology. This structure is 
highly suggestive of a 'Mach disk' commo nly seen in bow 
shock s associated with Herbig-Haro jets (jHartigan et al.l 
I1999H . In the flattened region, we find N e « 1600 cm" 3 , 
a factor «3 times greater than the environment. This is 
less than the maximum compression ratio of 4 expected 
for a strong shock, but is compatible with the Mach num- 
ber inferred from the shock speed (see below). 

From the high density region, we can determine a shock 
speed using the multitude of line ratios extracted from 
a single IFS exposure. Fig. l(m) shows theoretical line 
ratios from the radiative shock model oflHartiga n et all 
(|1987f ). The dashed horizontal lines indicate the maxi- 
mum and minimum shock speeds obtained from the new 
data for each corresponding line ratio. From the pres- 
ence of [Om], we infer a lower limit v s > 80 kms. 
[Ann] /He i suggests v s < 190 kms -1 , while both the 
[Sn]/Ha and [Arm] /Ha line ratios further constrain the 
shock velocity to v s ss 90 — 100 kms - . Thus, we esti- 
mate a shock speed v s ks 90 kms -1 , as indicated by the 
black dashed vertical line in Fig. l(m). This is broadly 
consistent with the line profiles being marginally resolved 
at w 120 - 150 km s -1 FWHM within the high density 
region. It also implies a Mach number of « 6, assum- 
ing an isothermal sound speed ~ 14 kms -1 for an ambi- 
ent gas temperature of 10 4 K. We note, however, that 
the deduced valu e of v s is likely to be somewhat over- 
estimated, as the lHartigan et all (|1987l ) shock model as- 
sumes a neutral pre-shock medium, and this produces a 
spectrum similar to preionized, shocked material with a 
lower v s (|Cox fc Ravmondl H985). 

3.3. Energy Budget 

The streamers are unlike the uniform wind-blown 
bow shocks associated with OB-runaway stars (see 
e.g. iKaper et alj Il997t ). Indeed, the the inferred 
space velocity of the 07 m companion is too low 
to be consistent with this interpretation. We esti- 
mate v* ~ 0.4 (M W /1Q- 7 M yr -1 ) 1 / 2 ^ 2 500 (p a /6 x 
10- 22 gcm- 3 ) -1 / 2 i?J.2kms -1 (IKaper et al.l fl99"l . 

where v w = lSOOu^isoo and M w are the wind velocity 
and mass loss rates of star 32 (|Pakull fc Angebaulti 
1986), p a is the mass density of the ambient medium, 
deduced from the [Sn] lines (§3.2), using an ioniza- 
tion fraction of x w 1, and where R = 3.2 R.3.2 is 
the distance between star 32 and the stagnation point 
(apex), which has an angular distance 12" corresponding 
to > 3.2 pc, where the inequality takes into account 
projection effects. We henceforth consider the possibility 

that the shock is produced by a jet from LMCX -1. 

Adopting the "dentist drill" model of lScheuerj (|1974f ). 
a cavity, or cocoon, develops naturally around a jet 
ploughing through the ambient ISM. A shock develops 
around the cocoon and the shocked gas then comes 
into dynamical (ram) pressure equilibrium with the 



jet. Thus, pjVj Ri p s Vg , where pj and p s are the mass 
densities of the jet and shocked gas, respectively, and 
Vj and v s are the corresponding speeds. The total 
jet kinetic power is Pj w n((j)jZj) 2 pjV^ (assuming it is 
not too relativistic), where fa is the jet half opening 
angle and Zj is the jet length (distance to the im- 
pact site, at the apex of the cocoon, from LMCX-1). 

Thus, Pj ~ 2 x 10 39 Vj i0 .i (fa, 10 ^,3.2 w s ,9o) 2 (p«/2 x 
10 -21 gem -3 ) ergs -1 . Here, fa.w = fa/ 10° 
([Miller- Jones. Fender fc Nakari 12009 find an up- 
per limit of fa < 10° for X-ray binary jets), 
z j,3.2 = z j sin -1 8/(3.2 pc) and 9 is the jet inclination 
angle to our line-of-sight. We have used a post-shock 
electron number density N e w 1.6 x 10 3 cm -3 (§3.2) to 
estimate p s . This estimate of jet power falls squarely 
within the range Pj w 10 37-42 erg s -1 estimated by 
iFender fc Poolevl (2000) for quasi-steady X-ray binary 
jets. 

It is noteworthy, however, that LMCX-1 currently ap- 
pears to be in a pers istent high/soft X-ray spectral state 
(|Schlegel et al.|[l994l ). du ring which it is unlikely to be 
producing a powerful jet (jFender. Belloni fc Galloll2"004h 
and its X-ray luminosity is dominated by the accretion 
disk, which radiates away all the accretion power. Since 
the accretion power inferred from the observed X-ray lu- 
minosity, « 2 x 10 38 ergs -1 , is an order of magnitude 
less than the inferred jet power, we sugges t that rather 
than being a dark jet (as lGallo et al.|[2005l proposed for 
CygX-1, which is in a low/hard spectral state), the non- 
detection of the LMCX-1 jet results from the inherent 
transient nature of X-ray binary jets. The jet is not 
presently seen because it may have recently switched 
off and the shock-ionized nebula surrounding LMCX- 
1 is still radiating the energy impacted during an ear- 
lier jet-active epoch. We estimate a radiative cooling 
timesc ale of just ~ 10 2 yrs, assu ming 0.3 solar metal- 
licity (D opita fc Sutherl and 2003]). The jet-active phase 
may have been triggered by a relatively sudden increase 
in the mass accretion rate, which has a present value 
M a ~ YlLx/c 2 w 4 x 10 -8 M Q yr -1 (assuming a ra- 
diative efficiency of for standard thin disk accre- 
tion). The dynamical timescale over which the jet-driven 
shock has been expanding into the ISM is estimated as 
tdyn ~ z j/v s ^ 0.03 Myr, which is remarkably similar to 
j et lifetime estimates obtained for o ther X-ray binaries 
(|Gallo et al.ll2005t iKaiser et al.ll200l . 

4. CONCLUSIONS 

We have demonstrated the extraordinary power of in- 
tegral field spectroscopy in our study of the emission-line 
nebula N159F around LMCX-1. Different emission-line 
diagnostics have been combined to separate out differ- 
ent ionization regions on extended scales. We propose 
that three ionizing sources contribute to the complex 
nature of N159F: 1. photoionization by star 32 (the 
07 III companion to LMCX-1); 2. X-ray photoioniza- 
tion by LMCX-1 itself; and 3. jet-driven shock ioniza- 
tion. The first two phenomena are relatively localized, 
whereas shock ionization dominates over the extended 
streamers, which appear to be bounded by a bow-shock 
morphology The morphology and ionization properties of 
the extended nebula around LMCX-1 are strikingly sim- 
ilar to those seen in other known jet-driven shock ionized 
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nebulae around Galactic X-ray binaries. We infer a shock 
speed of w s «90 km s _1 and deduce the power of the jet 
driving the shock to be w 2 x 10 39 erg s -1 . Since LMC X- 
1 is unlikely to be producing a jet in its current X-ray 
spectral state, we propose that the jet is intermittent, 
with a characteristic lifetime of ~ 0.03 Myr. 
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Fig. 1. — (a)-(d) Continuum subtracted integrated line flux images of N159F in He*, [Arm], [Ol], and [Om]; filled blue circles indi- 
cate the locations of stars 32 (right) and R148 (left); north is up and east is to the left (5" ~ 1.3 pc); (e) A diagnostic scatterplot of 
[O in]A5007/[O 1] A6300 versus [Ar m]A7136/[S n]A(6716+6731); the colour coding corresponds to different ionization regions shown in (f); 
(g) An ele ctron density ima ge, showing an arc of higher electron density to the south; (h) The Spitzer three colour composite image 
of N159F l| Jones et al.ll2005l ). constructed from the IRAC 3.6, 4.5 and 8.0 /im bands; (i)-(j) Ha and [Om]/[0 1] images overlaid with a 
proposed trail of the streamer filaments; (k) A false-colour image (R, G, B = [N 11] A6583, Ha, [O I1A630 0) with the proposed bow shock and 
jet orientation indicated; (1) Same as (f), with overlaid blue lines indicating the echelle slit positions of Ramsey ct al. (2006); grey lines arc 
the proposed streamers; (m) Predicted line ratios and corresponding shock speeds from the radiative shock model in Hartigan et al. (1987) 
(solid curves) overplotted with observed ranges of line ratios (dashed lines) measured from the region of high electron density shown in (g). 



